Legionnaires’ Disease Outbreaks
Holly Kessler
August 10, 2016
The coordinates of the Legionnairess Disease outbreaks will be used to gather weather data from surrounding stations. The averages of the data will be taken and outputted into a graph containing data from the last 10 years before the outbreak.
First I installed and downloaded the packages needed including rnoaa from github.
library(devtools)
install_github("geanders/rnoaa")
install_github("ropenscilabs/riem")library(devtools)
library(rnoaa)
library(riem)
library(countyweather)
library(dplyr)
library(plyr)
library(tidyr)
library(weathermetrics)
library(ggplot2)
library(lubridate)I created a data frame including the locations of each outbreak. I found the coordinates at http://maps.cga.harvard.edu/gpf/ and crossed checked them with Google coordinates. The other data in this set are year of outbreak and the year 10 years before the outbreak, onset date, and 14 days before the onset date.
| id | file_id | latitude | longitude | year_min | date_min | year_max | date_max | onset | before_onset |
|---|---|---|---|---|---|---|---|---|---|
| portugal | portugal | 38.96 | -8.99 | 2004 | 2004-01-01 | 2014 | 2014-12-31 | 2004-10-14 | 2004-09-30 |
| pittsburgh | pittsburgh | 40.43 | -79.98 | 2002 | 2002-01-01 | 2012 | 2012-12-31 | 2012-08-26 | 2012-08-12 |
| quebec | quebec | 46.85 | -71.34 | 2002 | 2002-01-01 | 2012 | 2012-12-31 | 2012-07-18 | 2012-07-04 |
| stoke-on-trent | stoke_on_trent | 53.02 | -2.15 | 2002 | 2002-01-01 | 2012 | 2012-12-31 | 2012-07-02 | 2012-06-18 |
| edinburgh | edinburgh | 55.94 | -3.20 | 2002 | 2002-01-01 | 2012 | 2012-12-31 | 2012-05-01 | 2012-04-17 |
| miyazaki | miyazaki | 31.89 | 131.34 | 1992 | 1992-01-01 | 2002 | 2002-12-31 | 2002-07-18 | 2002-07-04 |
| pas-de-calais | pas_de_calais | 50.51 | 2.37 | 1993 | 1993-01-01 | 2003 | 2003-12-31 | 2003-11-28 | 2003-11-14 |
| pamplona | pamplona | 42.81 | -1.65 | 1996 | 1996-01-01 | 2006 | 2006-12-31 | 2006-06-01 | 2006-05-18 |
| rapid city | rapid_city | 44.06 | -103.22 | 1995 | 1995-01-01 | 2005 | 2005-12-31 | 2005-05-26 | 2005-05-12 |
| sarpsborg | sarpsborg | 59.28 | 11.08 | 1995 | 1995-01-01 | 2005 | 2005-12-31 | 2005-05-12 | 2005-04-28 |
| barrow-in-furness | barrow_in_furness | 54.10 | -3.22 | 1992 | 1992-01-01 | 2002 | 2002-12-31 | 2002-07-30 | 2002-07-16 |
| murcia | murcia | 37.98 | -1.12 | 1991 | 1991-01-01 | 2001 | 2001-12-31 | 2001-06-26 | 2001-06-12 |
| melbourne | melbourne | -37.86 | 145.07 | 1990 | 1990-01-01 | 2000 | 2000-12-31 | 2000-04-17 | 2000-04-03 |
| bovenkarspel | bovenkarspel | 52.70 | 5.24 | 1989 | 1989-01-01 | 1999 | 1999-12-31 | 1999-02-25 | 1999-02-11 |
| london | london | 51.52 | -0.10 | 1979 | 1979-01-01 | 1989 | 1989-12-31 | 1989-01-01 | 1988-12-18 |
| sydney | sydney | -33.85 | 150.93 | 2006 | 2006-01-01 | 2016 | 2016-12-31 | 2016-04-25 | 2016-04-11 |
| genesee1 | genesee1 | 43.09 | -83.63 | 2004 | 2004-01-01 | 2014 | 2014-12-31 | 2014-06-06 | 2014-05-23 |
| genesee2 | genesee2 | 43.09 | -83.63 | 2005 | 2005-01-01 | 2015 | 2015-12-31 | 2015-05-04 | 2015-04-20 |
| columbus | columbus | 39.98 | -82.99 | 2003 | 2003-01-01 | 2013 | 2013-12-31 | 2013-07-09 | 2013-06-25 |
| bronx | bronx | 40.82 | -73.92 | 2005 | 2005-01-01 | 2015 | 2015-12-31 | 2015-07-12 | 2015-06-28 |
The next function will download information from all of the stations. It only needs to be downloaded once per session. It will take a couple minutes to download.
I created a loop to get a list of the stations within 30 km of the location.
station_data <- ghcnd_stations()[[1]]
df <- list()
for(i in 1:length(outbreak_loc$id))
{
df[[i]] <- (meteo_nearby_stations(lat_lon_df = outbreak_loc[i,],
station_data = station_data,
var = c("PRCP","TAVG","TMAX","TMIN",
"AWND","MDPR"),
year_min = outbreak_loc[i, "year_min"],
year_max = outbreak_loc[i, "year_max"],
radius = 30))
}
names(df) <- outbreak_loc$id
stations <- df
saveRDS(stations, file = "stations.RData")stations <- readRDS("stations.RData")
stations## $`barrow-in-furness`
## [1] id name latitude longitude distance
## <0 rows> (or 0-length row.names)
##
## $bovenkarspel
## id name latitude longitude distance
## 1 NLE00101917 ENKHUIZEN 52.6917 5.2944 3.780363
## 2 NLE00109144 HOOGKARSPEL 52.6867 5.1669 5.143626
## 3 NLE00101928 MEDEMBLIK 52.7781 5.1014 12.746887
## 4 NLE00100501 HOORN 52.6444 5.0681 13.136277
## 5 NLE00102479 BERKHOUT 52.6428 4.9789 18.718869
## 6 NLE00109146 HOOGWOUD 52.7281 4.9608 19.065006
## 7 NLE00109174 KREILEROORD 52.8619 5.0953 20.464698
## 8 NLE00102134 STAVOREN 52.8967 5.3831 23.894441
## 9 NLE00109232 OBDAM 52.6775 4.8769 24.600539
## 10 NLE00109054 EDAM 52.5114 5.0467 24.701897
## 11 NLE00109250 OUDEMIRDUM 52.8608 5.5078 25.379621
## 12 NLE00109162 KOLHORN 52.7914 4.8919 25.540549
## 13 NLE00109354 WEST BEEMSTER 52.5817 4.9028 26.281190
## 14 NLE00101948 TOLLEBEEK 52.6719 5.6300 26.472774
## 15 NLE00101930 DEN OEVER 52.9217 5.0383 28.133541
## 16 NLE00101932 MARKEN 52.4600 5.1078 28.142012
## 17 NLE00109028 DE HAUKES 52.8783 4.9408 28.246661
## 18 NLE00109254 PURMEREND 52.5125 4.9506 28.575998
##
## $bronx
## id name latitude longitude distance
## 1 USW00014732 NY NEW YORK LAGUARDIA AP 40.7794 -73.8803 5.616756
## 2 USC00305799 NY NEW YORK BOTANICAL GRD 40.8667 -73.8833 6.041176
## 3 USW00094728 NY NEW YORK CNTRL PK TWR 40.7789 -73.9692 6.167420
## 4 USC00300961 NY BRONX 40.8369 -73.8494 6.230297
## 5 US1NJBG0018 NJ PALISADES PARK 0.2 WNW 40.8481 -74.0002 7.435652
## 6 USC00305804 NY NEW YORK LAUREL HILL 40.7333 -73.9333 9.705430
## 7 US1NJBG0003 NJ TENAFLY 1.3 W 40.9147 -73.9775 11.587163
## 8 US1NYQN0002 NY MIDDLE VILLAGE 0.5 SW 40.7145 -73.8819 12.161951
## 9 USW00094741 NJ TETERBORO AP 40.8500 -74.0614 12.354793
## 10 US1NJBG0001 NJ BERGENFIELD 0.3 SW 40.9213 -74.0020 13.206762
## 11 US1NJBG0012 NJ WOOD RIDGE 0.6 SE 40.8420 -74.0830 13.930427
## 12 USC00284339 NJ JERSEY CITY 40.7419 -74.0572 14.452190
## 13 USC00284931 NJ LODI 40.8667 -74.0833 14.685375
## 14 US1NJBG0033 NJ WOOD RIDGE 0.6 NNW 40.8536 -74.0943 15.131881
## 15 US1NYWC0009 NY NEW ROCHELLE 1.3 S 40.9040 -73.7770 15.226895
## 16 US1NJBG0013 NJ RUTHERFORD 1.2 N 40.8373 -74.1065 15.809146
## 17 US1NYKN0025 NY BROOKLYN 3.1 NW 40.6846 -73.9867 16.069963
## 18 US1NJBG0031 NJ DEMAREST 0.6 NNW 40.9628 -73.9600 16.230719
## 19 US1NJBG0002 NJ SADDLE BROOK TWP 0.6 E 40.9027 -74.0834 16.534408
## 20 US1NJBG0011 NJ NORTH ARLINGTON 0.7 NE 40.7944 -74.1190 16.988985
## 21 US1NJBG0008 NJ SADDLE BROOK TWP 0.3 NNE 40.9071 -74.0934 17.505118
## 22 USC00286146 NJ NEW MILFORD 40.9611 -74.0158 17.635536
## 23 US1NJBG0015 NJ NORTH ARLINGTON 0.7 WNW 40.7915 -74.1398 18.769309
## 24 US1NJHD0002 NJ KEARNY 1.7 NW 40.7729 -74.1409 19.318492
## 25 US1NJBG0005 NJ WESTWOOD 0.8 ESE 40.9830 -74.0159 19.836072
## 26 US1NJBG0010 NJ RIVER VALE TWP 1.5 S 40.9915 -74.0123 20.587159
## 27 US1NYNS0007 NY FLORAL PARK 0.4 W 40.7230 -73.7110 20.642010
## 28 USC00307497 NY SCARSDALE 40.9833 -73.8000 20.770959
## 29 US1NJHD0001 NJ HARRISON 0.3 N 40.7480 -74.1518 21.094543
## 30 USC00286775 NJ PATERSON 40.9000 -74.1500 21.289996
## 31 USC00283704 NJ HARRISON 40.7514 -74.1567 21.338273
## 32 US1NJBG0020 NJ PARAMUS 1.8 NNW 40.9682 -74.0902 21.822558
## 33 USC00302129 NY DOBBS FERRY-ARDSLEY 41.0072 -73.8344 22.023427
## 34 US1NJBG0017 NJ GLEN ROCK 0.7 SSE 40.9511 -74.1183 22.145027
## 35 US1NJES0020 NJ BLOOMFIELD 1.7 S 40.7850 -74.1885 22.932509
## 36 US1NYKN0003 NY BROOKLYN 2.4 SW 40.6194 -73.9859 22.986706
## 37 US1NYWC0005 NY HARRISON 4.1 SSW 40.9639 -73.7232 23.014851
## 38 USC00307587 NY SEA CLIFF 40.8506 -73.6483 23.109762
## 39 US1NJBG0037 NJ GLEN ROCK 0.4 WNW 40.9614 -74.1328 23.815579
## 40 US1NJPS0014 NJ HAWTHORNE 1.0 SSE 40.9436 -74.1523 23.880745
## 41 USC00289832 NJ WOODCLIFF LAKE 41.0139 -74.0425 23.891667
## 42 US1NJES0015 NJ MONTCLAIR 2.2 NNE 40.8565 -74.2004 23.935385
## 43 USW00094789 NY NEW YORK JFK INTL AP 40.6386 -73.7622 24.159135
## 44 US1NJPS0017 NJ WOODLAND PARK 0.1 NW 40.8918 -74.1960 24.547061
## 45 USC00305382 NY MINEOLA 1 W 40.7333 -73.6500 24.694525
## 46 US1NJPS0005 NJ HAWTHORNE 0.4 S 40.9519 -74.1577 24.787058
## 47 US1NJES0011 NJ CEDAR GROVE TWP 0.9 NE 40.8648 -74.2157 25.368252
## 48 USC00305796 NY NY AVE V BROOKLYN 40.5939 -73.9808 25.658206
## 49 US1NJPS0018 NJ PATERSON 2.0 W 40.9163 -74.2005 25.903423
## 50 USW00014734 NJ NEWARK INTL AP 40.6825 -74.1694 25.982975
## 51 US1NJPS0003 NJ LITTLE FALLS TWP 0.2 NNW 40.8788 -74.2205 26.107408
## 52 USC00305380 NY MINEOLA 1 NE 40.7494 -73.6233 26.184670
## 53 USC00303786 NY HEMPSTEAD MALVERNE 40.6833 -73.6667 26.197491
## 54 USW00014786 NY NEW YORK FLOYD BENNETT FLD 40.5833 -73.8833 26.501042
## 55 USC00305821 NY NY WESTERLEIGH STAT IS 40.6333 -74.1167 26.565419
## 56 US1NJPS0012 NJ LITTLE FALLS TWP 0.5 WNW 40.8796 -74.2270 26.658883
## 57 USC00303781 NY HEMPSTEAD GARDEN CITY 40.7167 -73.6333 26.737225
## 58 USC00284887 NJ LITTLE FALLS 40.8858 -74.2261 26.764590
## 59 US1NJES0024 NJ CEDAR GROVE TWP 0.4 W 40.8557 -74.2356 26.845270
## 60 US1NYNS0014 NY LYNBROOK 0.3 NW 40.6623 -73.6780 26.891776
## 61 USC00285503 NJ MIDLAND PARK 40.9939 -74.1453 27.062890
## 62 USC00305377 NY MINEOLA 40.7328 -73.6183 27.191818
## 63 US1NJPS0004 NJ NORTH HALEDON 0.6 N 40.9713 -74.1856 27.953822
## 64 US1NJES0010 NJ VERONA TWP 0.7 SW 40.8255 -74.2531 28.035405
## 65 US1NJES0021 NJ VERONA TWP 0.6 WSW 40.8305 -74.2539 28.119240
## 66 US1NJES0004 NJ NORTH CALDWELL 0.6 SSE 40.8576 -74.2523 28.265572
## 67 USC00309117 NY WESTBURY 40.7333 -73.6000 28.617850
## 68 US1NJPS0008 NJ WAYNE TWP 1.1 ESE 40.9412 -74.2267 29.094312
## 69 USC00309270 NY WEST NYACK 41.0833 -73.9667 29.539139
## 70 US1NYWC0003 NY WHITE PLAINS 3.1 NNW 41.0639 -73.7722 29.826697
## 71 US1NYNS0009 NY MILL NECK 1.1 SW 40.8704 -73.5717 29.828998
## 72 US1NYRL0005 NY WEST NYACK 1.3 WSW 41.0835 -73.9930 29.934369
##
## $columbus
## id name latitude longitude
## 1 USW00093891 OH COLUMBUS WSO CITY 39.9667 -83.0000
## 2 US1OHFR0018 OH COLUMBUS 2.4 WNW 39.9977 -83.0323
## 3 US1OHFR0025 OH COLUMBUS 2.8 WSW 39.9804 -83.0397
## 4 US1OHFR0003 OH GRANDVIEW HEIGHTS 0.1 N 39.9810 -83.0401
## 5 USC00331782 OH COLUMBUS UNIV FARM 40.0000 -83.0500
## 6 USC00331785 OH COLUMBUS WCMH 40.0250 -83.0269
## 7 US1OHFR0039 OH COLUMBUS 3.4 WNW 40.0143 -83.0434
## 8 US1OHFR0034 OH COLUMBUS 3.6 NW 40.0191 -83.0437
## 9 US1OHFR0020 OH COLUMBUS 3.5 NE 40.0287 -82.9477
## 10 US1OHFR0001 OH UPPER ARLINGTON 0.9 E 40.0279 -83.0543
## 11 US1OHFR0021 OH MARBLE CLIFF 1.1 WNW 39.9931 -83.0786
## 12 US1OHFR0007 OH UPPER ARLINGTON 1.3 SSW 40.0112 -83.0832
## 13 USW00014821 OH COLUMBUS PORT COLUMBUS INTL AP 39.9914 -82.8808
## 14 USC00331783 OH COLUMBUS-VALLEY CROSSING 39.9047 -82.9200
## 15 USC00331781 OH COLUMBUS SULLIVANT AVE 39.9500 -83.1167
## 16 US1OHFR0012 OH UPPER ARLINGTON 2.4 NNW 40.0604 -83.0815
## 17 USC00331777 OH COLUMBUS-HAP CREMEAN WP 40.0603 -82.8942
## 18 US1OHFR0024 OH COLUMBUS 9.3 NNE 40.0925 -82.9582
## 19 USW00004804 OH COLUMBUS OHIO STATE UNIV AP 40.0781 -83.0781
## 20 US1OHFR0037 OH REYNOLDSBURG 1.6 W 39.9588 -82.8294
## 21 US1OHFR0016 OH DUBLIN 3.7 ESE 40.0923 -83.0725
## 22 US1OHFR0022 OH GALLOWAY 3.1 N 39.9561 -83.1592
## 23 USC00331779 OH COLUMBUS-PARSONS AVE. 39.8469 -82.9872
## 24 USC00338951 OH WESTERVILLE 40.1264 -82.9433
## 25 US1OHFR0010 OH WESTERVILLE 0.2 WNW 40.1226 -82.9213
## 26 US1OHFR0030 OH HILLIARD 1.8 W 40.0344 -83.1768
## 27 US1OHFR0008 OH NEW ALBANY 2.8 SSE 40.0403 -82.7980
## 28 US1OHFR0002 OH DUBLIN 3.2 ENE 40.1299 -83.0742
## 29 USW00013812 OH COLUMBUS RICKENBACKER 39.8167 -82.9333
## 30 US1OHLC0002 OH PATASKALA 4.4 WNW 40.0273 -82.7490
## 31 US1OHFF0010 OH PICKERINGTON 0.1 SSE 39.8895 -82.7680
## 32 US1OHFF0005 OH PICKERINGTON 2.7 NNE 39.9263 -82.7469
## 33 US1OHDL0002 OH WESTERVILLE 4.0 N 40.1790 -82.9256
## 34 US1OHDL0012 OH LEWIS CENTER 1.4 E 40.1967 -82.9835
## 35 US1OHDL0008 OH DUBLIN 3.7 NNW 40.1642 -83.1687
## 36 US1OHFR0023 OH HARRISBURG 3.7 WNW 39.8378 -83.2321
## 37 USC00334616 OH LITHOPOLIS 2 S 39.7833 -82.8167
## 38 USC00330083 OH ALEXANDRIA 4 WSW 40.0833 -82.6833
## 39 US1OHDL0013 OH GALENA 3.0 WNW 40.2337 -82.9362
## 40 US1OHLC0011 OH PATASKALA 2.0 NE 40.0240 -82.6511
## 41 USC00332170 OH DERBY 39.7667 -83.2000
## distance
## 1 1.706827
## 2 4.106137
## 3 4.234920
## 4 4.270197
## 5 5.574388
## 6 5.909011
## 7 5.936178
## 8 6.310790
## 9 6.504242
## 10 7.639623
## 11 7.687720
## 12 8.664212
## 13 9.389591
## 14 10.282003
## 15 11.301408
## 16 11.858973
## 17 12.094496
## 18 12.799033
## 19 13.238469
## 20 13.887609
## 21 14.326951
## 22 14.662096
## 23 14.801971
## 24 16.757176
## 25 16.900247
## 26 17.021110
## 27 17.673463
## 28 18.143413
## 29 18.791310
## 30 21.190424
## 31 21.436847
## 32 21.564691
## 33 22.796070
## 34 24.102284
## 35 25.509352
## 36 26.008132
## 37 26.401695
## 38 28.527330
## 39 28.578803
## 40 29.278406
## 41 29.727063
##
## $edinburgh
## [1] id name latitude longitude distance
## <0 rows> (or 0-length row.names)
##
## $genesee1
## id name latitude longitude distance
## 1 USC00201150 MI BURTON 4N 43.0675 -83.5919 3.979309
## 2 US1MIGN0010 MI BURTON 0.9 NNW 43.0085 -83.6274 9.064849
## 3 US1MIGN0008 MI MOUNT MORRIS 3.1 WSW 43.1057 -83.7580 10.538333
## 4 US1MIGN0014 MI DAVISON 3.3 SW 43.0003 -83.5684 11.159857
## 5 US1MIGN0005 MI DAVISON 0.7 SSW 43.0219 -83.5246 11.431375
## 6 USC00202851 MI FLINT 7 W 43.0378 -83.7694 12.725462
## 7 US1MIGN0024 MI CLIO 0.4 SW 43.1725 -83.7423 12.930649
## 8 USC00201645 MI CLIO 43.1794 -83.7369 13.193328
## 9 US1MIGN0009 MI BURTON 3.3 SW 42.9613 -83.6636 14.569099
## 10 USW00014826 MI FLINT BISHOP INTL AP 42.9667 -83.7494 16.797891
## 11 US1MIGN0023 MI GRAND BLANC 3.8 WNW 42.9440 -83.6886 16.919078
## 12 US1MISG0004 MI BIRCH RUN 2.6 ESE 43.2291 -83.7470 18.146491
## 13 US1MIGN0020 MI SWARTZ CREEK 2.0 NNE 42.9897 -83.8166 18.824526
## 14 US1MIGN0015 MI FLINT 6.4 SSW 42.9326 -83.7210 19.001794
## 15 US1MIGN0018 MI GRAND BLANC 0.7 SE 42.9187 -83.6079 19.132279
## 16 USC00203278 MI GOODRICH 42.9164 -83.5097 21.640755
## 17 USC00204659 MI LAPEER 2W 43.0581 -83.3606 22.167570
## 18 US1MIGN0022 MI GRAND BLANC 2.9 SE 42.8909 -83.5858 22.428899
## 19 USC00204662 MI LAPEER 4 SW 43.0081 -83.3750 22.633695
## 20 US1MIGN0004 MI MONTROSE 0.4 NW 43.1794 -83.8987 23.962697
## 21 USC00205488 MI MILLINGTON 3 SE 43.2836 -83.4792 24.756898
## 22 US1MILP0003 MI LAPEER 1.1 SSW 43.0316 -83.3293 25.277892
## 23 USC00201721 MI COLUMBIAVILLE 5 NE 43.1897 -83.3358 26.319366
## 24 USC00204655 MI LAPEER WWTP 43.0608 -83.3075 26.394851
## 25 USC00202955 MI FRANKENMUTH 1SE 43.3194 -83.7161 26.445482
##
## $genesee2
## id name latitude longitude distance
## 1 USC00201150 MI BURTON 4N 43.0675 -83.5919 3.979309
## 2 US1MIGN0010 MI BURTON 0.9 NNW 43.0085 -83.6274 9.064849
## 3 US1MIGN0008 MI MOUNT MORRIS 3.1 WSW 43.1057 -83.7580 10.538333
## 4 US1MIGN0014 MI DAVISON 3.3 SW 43.0003 -83.5684 11.159857
## 5 US1MIGN0005 MI DAVISON 0.7 SSW 43.0219 -83.5246 11.431375
## 6 USC00202851 MI FLINT 7 W 43.0378 -83.7694 12.725462
## 7 US1MIGN0024 MI CLIO 0.4 SW 43.1725 -83.7423 12.930649
## 8 USC00201645 MI CLIO 43.1794 -83.7369 13.193328
## 9 US1MIGN0009 MI BURTON 3.3 SW 42.9613 -83.6636 14.569099
## 10 USW00014826 MI FLINT BISHOP INTL AP 42.9667 -83.7494 16.797891
## 11 US1MIGN0023 MI GRAND BLANC 3.8 WNW 42.9440 -83.6886 16.919078
## 12 US1MISG0004 MI BIRCH RUN 2.6 ESE 43.2291 -83.7470 18.146491
## 13 US1MIGN0020 MI SWARTZ CREEK 2.0 NNE 42.9897 -83.8166 18.824526
## 14 US1MIGN0015 MI FLINT 6.4 SSW 42.9326 -83.7210 19.001794
## 15 US1MIGN0018 MI GRAND BLANC 0.7 SE 42.9187 -83.6079 19.132279
## 16 USC00203278 MI GOODRICH 42.9164 -83.5097 21.640755
## 17 USC00204659 MI LAPEER 2W 43.0581 -83.3606 22.167570
## 18 US1MIGN0022 MI GRAND BLANC 2.9 SE 42.8909 -83.5858 22.428899
## 19 USC00204662 MI LAPEER 4 SW 43.0081 -83.3750 22.633695
## 20 US1MIGN0004 MI MONTROSE 0.4 NW 43.1794 -83.8987 23.962697
## 21 USC00205488 MI MILLINGTON 3 SE 43.2836 -83.4792 24.756898
## 22 US1MILP0003 MI LAPEER 1.1 SSW 43.0316 -83.3293 25.277892
## 23 USC00201721 MI COLUMBIAVILLE 5 NE 43.1897 -83.3358 26.319366
## 24 USC00204655 MI LAPEER WWTP 43.0608 -83.3075 26.394851
## 25 USC00202955 MI FRANKENMUTH 1SE 43.3194 -83.7161 26.445482
##
## $london
## id name latitude longitude distance
## 1 UKM00003772 HEATHROW 51.478 -0.461 25.42177
##
## $melbourne
## id name latitude longitude
## 1 ASN00086141 BURWOOD PLC -37.8333 145.0833
## 2 ASN00086018 CAULFIELD (RACECOURSE) -37.8795 145.0368
## 3 ASN00086304 HAWTHORN (SCOTCH COLLEGE) -37.8361 145.0294
## 4 ASN00086019 CAULFIELD SOUTH -37.9000 145.0333
## 5 ASN00086060 KEW -37.8167 145.0333
## 6 ASN00086095 PRAHRAN (COMO HOUSE) -37.8376 145.0048
## 7 ASN00086088 OAKLEIGH (METROPOLITAN GOLF CL -37.9142 145.0935
## 8 ASN00086012 BOX HILL AGED MENS RETREAT -37.8364 145.1364
## 9 ASN00086002 BALWYN NORTH -37.8000 145.0500
## 10 ASN00086016 CANTERBURY -37.8000 145.0500
## 11 ASN00086222 SYNDAL CSIRO -37.8833 145.1500
## 12 ASN00086006 BENTLEIGH -37.9279 145.0749
## 13 ASN00086011 BOX HILL SOUTH -37.8333 145.1500
## 14 ASN00086033 BRIGHTON (DENDY PARK BOWLING C -37.9252 145.0254
## 15 ASN00086340 BLACKBURN -37.8167 145.1500
## 16 ASN00086232 MELBOURNE BOTANICAL GARDENS -37.8303 144.9767
## 17 ASN00086013 BURNLEY -37.8000 145.0000
## 18 ASN00086100 BURNLEY RICHMOND PARK -37.8000 145.0000
## 19 ASN00086307 ABBOTSFORD -37.8000 145.0000
## 20 ASN00086341 DONCASTER -37.7833 145.1167
## 21 ASN00086279 NORTHCOTE -37.7797 145.0314
## 22 ASN00086316 VERMONT TRANSPORT RESEARCH -37.8587 145.1847
## 23 ASN00086071 MELBOURNE REGIONAL OFFICE -37.8075 144.9700
## 24 ASN00086020 CHELTENHAM KINGSTON CENTRE -37.9551 145.0782
## 25 ASN00086317 SANDRINGHAM 2 -37.9500 145.0167
## 26 ASN00086031 DONCASTER -37.7833 145.1500
## 27 ASN00086303 GLEN WAVERLEY (GOLF COURSE) -37.8886 145.1928
## 28 ASN00086074 MITCHAM -37.8219 145.1906
## 29 ASN00086257 VERMONT SOUTH -37.8500 145.2000
## 30 ASN00086260 HEIDELBERG MMBW -37.7567 145.0533
## 31 ASN00086348 MULGRAVE -37.9167 145.1833
## 32 ASN00086047 GLEN WAVERLEY (SYNDAL) -37.9000 145.2000
## 33 ASN00086241 TEMPLESTOWE PASSIONATE FATHERS -37.7667 145.1500
## 34 ASN00086302 WANTIRNA SOUTH -37.8667 145.2167
## 35 ASN00086378 BRUNSWICK -37.7667 144.9797
## 36 ASN00086111 SPRINGVALE NECROPOLIS -37.9445 145.1764
## 37 ASN00086369 SPRINGVALE (SANDOWN) -37.9535 145.1655
## 38 ASN00086068 VIEWBANK (ARPANSA) -37.7408 145.0972
## 39 ASN00086077 MOORABBIN AIRPORT -37.9800 145.0964
## 40 ASN00086146 BEAUMARIS -37.9771 145.0273
## 41 ASN00086053 HEIDELBERG (ROSANNA) -37.7500 145.0000
## 42 ASN00086333 THORNBURY -37.7500 145.0000
## 43 ASN00087145 WILLIAMSTOWN BOWLING CLUB -37.8500 144.9000
## 44 ASN00086362 DONCASTER (MANNINGHAM DEPOT) -37.7494 145.1703
## 45 ASN00086024 COBURG -37.7500 144.9667
## 46 ASN00086263 MORDIALLOC NORTH -38.0000 145.1000
## 47 ASN00086351 BUNDOORA (LATROBE UNIVERSITY) -37.7163 145.0453
## 48 ASN00086039 FLEMINGTON RACECOURSE -37.7915 144.9067
## 49 ASN00086104 SCORESBY RESEARCH INSTITUTE -37.8710 145.2561
## 50 ASN00086096 PRESTON RESERVOIR -37.7214 145.0059
## 51 ASN00086230 BAYSWATER -37.8372 145.2558
## 52 ASN00086379 RINGWOOD NORTH -37.7917 145.2433
## 53 ASN00087022 FOOTSCRAY -37.8167 144.8833
## 54 ASN00086101 RINGWOOD -37.8000 145.2500
## 55 ASN00086313 WARRANDYTE -37.7469 145.2098
## 56 ASN00086284 LYSTERFIELD (ALBION REID) -37.9333 145.2500
## 57 ASN00086343 BUNDOORA -37.7000 145.0667
## 58 ASN00086122 WATSONIA LOYOLA -37.7000 145.0833
## 59 ASN00086347 YARRA RIVER @ WARRANDYTE -37.7417 145.2167
## 60 ASN00086091 PASCOE VALE -37.7333 144.9333
## 61 ASN00086224 DANDENONG -37.9785 145.2235
## 62 ASN00086035 ELTHAM -37.7011 145.1547
## 63 ASN00087038 MARIBYRNONG EXPLOSIVES FACTORY -37.7750 144.8767
## 64 ASN00086147 ASPENDALE CSIRO -38.0333 145.1000
## 65 ASN00086332 ESSENDON 2 -37.7333 144.9167
## 66 ASN00086046 GLENROY -37.7167 144.9333
## 67 ASN00086324 FERNTREE GULLY (PROBERT) -37.8797 145.2964
## 68 ASN00086027 CROYDON (SAMUEL STREET) -37.7903 145.2812
## 69 ASN00086010 BORONIA 1 -37.8500 145.3000
## 70 ASN00086215 DIAMOND CREEK -37.6833 145.1333
## 71 ASN00086234 CROYDON (COUNCIL DEPOT) -37.7869 145.2847
## 72 ASN00086038 ESSENDON AIRPORT -37.7276 144.9066
## 73 ASN00087131 ALTONA (CITY OFFICES) -37.8633 144.8261
## 74 ASN00086110 MORANG SOUTH -37.6667 145.0833
## 75 ASN00087148 SUNSHINE (CITY OF BRINBANK) -37.7928 144.8344
## 76 ASN00086355 PLENTY (DAPRANO) -37.6667 145.1333
## 77 ASN00086250 PLENTY -37.6600 145.1244
## 78 ASN00086286 GLENROY (EVELL STREET) -37.6833 144.9333
## 79 ASN00086256 FERNY CREEK -37.8833 145.3333
## 80 ASN00086292 BANGHOLME -38.0500 145.1833
## 81 ASN00086210 BONBEACH (CARRUM) -38.0651 145.1294
## 82 ASN00086372 FERNY CREEK (DUNNS HILL) -37.8775 145.3364
## 83 ASN00086228 KILSYTH (TUKAMUNDA) -37.8167 145.3333
## 84 ASN00086365 MOOROOLBARK -37.7792 145.3197
## 85 ASN00086251 UPWEY SHIRE COUNCIL -37.9144 145.3317
## 86 ASN00086270 SASSAFRAS (THE CRESCENT) -37.8667 145.3500
## 87 ASN00086221 SHERBROOKE BURNAM BEECHES -37.8833 145.3500
## 88 ASN00086243 MOUNT DANDENONG GTV9 -37.8306 145.3500
## 89 ASN00086059 KANGAROO GROUND -37.6830 145.2518
## 90 ASN00086258 BELGRAVE -37.9167 145.3500
## 91 ASN00086103 SASSAFRAS POST OFFICE -37.8000 145.3500
## 92 ASN00086254 CARRUM DOWNS SEWER WORKS -38.0783 145.1733
## 93 ASN00086036 EPPING -37.6312 144.9846
## 94 ASN00086287 HURSTBRIDGE (NORALA) -37.6500 145.2167
## 95 ASN00086066 LILYDALE -37.7488 145.3416
## 96 ASN00086076 MONTROSE -37.8019 145.3675
## 97 ASN00087020 DEER PARK -37.7667 144.7833
## 98 ASN00086085 NARRE WARREN NORTH (NARRE WARR -37.9897 145.3356
## 99 ASN00087031 LAVERTON RAAF -37.8565 144.7566
## 100 ASN00087027 KEILOR (ARUNDEL) -37.6942 144.8342
## 101 ASN00086056 KALLISTA FORESTRY -37.9000 145.3833
## 102 ASN00086086 NARRE WARREN (BRECHIN) -38.0333 145.3000
## 103 ASN00086319 DOREEN -37.6167 145.1500
## 104 ASN00087032 LAVERTON SALINES -37.9000 144.7500
## 105 ASN00087177 LAVERTON COMPARISON -37.8633 144.7456
## 106 ASN00086225 CARNEGIE -38.1167 145.0667
## 107 ASN00086305 GREENVALE RESERVOIR -37.6369 144.9072
## 108 ASN00086384 MELBOURNE AIRPORT COMPARISON -37.6750 144.8419
## 109 ASN00087015 KEILOR -37.7025 144.8072
## 110 ASN00087107 POINT COOK RAAF ACADEMY -37.9333 144.7500
## 111 ASN00087106 POINT COOK RAAF (DUNNINGS RD) -37.8833 144.7333
## 112 ASN00086107 SILVAN STATE SCHOOL -37.8000 145.4000
## distance
## 1 3.190327
## 2 3.632396
## 3 4.446429
## 4 5.491613
## 5 5.793822
## 6 6.243145
## 7 6.369850
## 8 6.393542
## 9 6.899047
## 10 6.899047
## 11 7.484784
## 12 7.562369
## 13 7.626099
## 14 8.238821
## 15 8.516796
## 16 8.833034
## 17 9.072309
## 18 9.072309
## 19 9.072309
## 20 9.463805
## 21 9.551015
## 22 10.070617
## 23 10.545355
## 24 10.599081
## 25 11.046224
## 26 11.050534
## 27 11.237859
## 28 11.406147
## 29 11.467486
## 30 11.579751
## 31 11.773203
## 32 12.245870
## 33 12.530662
## 34 12.899710
## 35 13.059611
## 36 13.245203
## 37 13.352636
## 38 13.468157
## 39 13.542852
## 40 13.548961
## 41 13.690472
## 42 13.690472
## 43 14.966642
## 44 15.129262
## 45 15.230624
## 46 15.788086
## 47 16.125457
## 48 16.239778
## 49 16.382188
## 50 16.408663
## 51 16.509685
## 52 17.010511
## 53 17.087517
## 54 17.158761
## 55 17.578887
## 56 17.773351
## 57 17.793552
## 58 17.829544
## 59 18.416453
## 60 18.513537
## 61 18.839479
## 62 19.172878
## 63 19.432886
## 64 19.448805
## 65 19.491447
## 66 19.954956
## 67 19.993316
## 68 20.103933
## 69 20.223613
## 70 20.420695
## 71 20.535018
## 72 20.564229
## 73 21.414594
## 74 21.525752
## 75 22.000521
## 76 22.202552
## 77 22.747356
## 78 23.030744
## 79 23.256162
## 80 23.345895
## 81 23.393049
## 82 23.465245
## 83 23.617812
## 84 23.701976
## 85 23.749365
## 86 24.591309
## 87 24.713433
## 88 24.802430
## 89 25.351527
## 90 25.367641
## 91 25.480101
## 92 25.907856
## 93 26.526359
## 94 26.675842
## 95 26.875070
## 96 26.914609
## 97 27.238327
## 98 27.399212
## 99 27.516720
## 100 27.737658
## 101 27.854473
## 102 27.894067
## 103 27.953385
## 104 28.435114
## 105 28.480727
## 106 28.545203
## 107 28.640881
## 108 28.725743
## 109 28.984935
## 110 29.237785
## 111 29.667485
## 112 29.740401
##
## $miyazaki
## id name latitude longitude distance
## 1 JA000047830 MIYAZAKI 31.933 131.417 8.699733
## 2 JA000047829 MIYAKONOJO 31.733 131.083 29.908209
##
## $murcia
## id name latitude longitude distance
## 1 SPE00120323 MURCIA 38.0028 -1.1692 5.001690
## 2 SPE00120332 MURCIA/ALCANTARILLA 37.9578 -1.2294 9.902609
##
## $pamplona
## id name latitude longitude distance
## 1 SPE00120350 PAMPLONA (OBSERVATORIO) 42.8175 -1.6364 1.387848
## 2 SPE00120359 PAMPLONA 42.7767 -1.6500 3.702791
##
## $`pas-de-calais`
## [1] id name latitude longitude distance
## <0 rows> (or 0-length row.names)
##
## $pittsburgh
## id name latitude longitude distance
## 1 US1PAAL0014 PA PITTSBURGH 1.6 SW 40.4226 -79.9974 1.687108
## 2 USW00014861 PA PITTSBURGH WSO CITY 40.4500 -80.0000 2.794735
## 3 USC00366986 PA PIT CARNEGIE SCI CTR 40.4450 -80.0186 3.667959
## 4 US1PAAL0017 PA WHITEHALL 1.0 SW 40.3475 -80.0022 9.364279
## 5 USW00014762 PA PITTSBURGH ALLEGHENY CO AP 40.3547 -79.9217 9.720301
## 6 US1PAAL0011 PA WEST MIFFLIN 1.3 SW 40.3466 -79.9283 10.255413
## 7 US1PAAL0031 PA SCOTT TOWNSHIP 1.3 NW 40.3978 -80.0967 10.508788
## 8 USC00360861 PA BRADDOCK LOCK 2 40.3917 -79.8594 11.063215
## 9 US1PAAL0009 PA PATHFINDER 40.3416 -80.0485 11.414110
## 10 USC00361773 PA CORAOPOLIS NEVILLE IS 40.5000 -80.0833 11.702719
## 11 USC00362574 PA EMSWORTH L/D OHIO RVR 40.5019 -80.0833 11.844197
## 12 US1PAAL0016 PA GLENSHAW 1.3 NW 40.5488 -79.9800 13.209957
## 13 USC00365573 PA MCKEESPORT 40.3392 -79.8603 14.308275
## 14 US1PAAL0008 PA UPPER ST. CLAIR 1.7 WNW 40.3412 -80.1026 14.329100
## 15 USC00361033 PA BRUCETON 1 S 40.3000 -79.9833 14.458044
## 16 US1PAAL0020 PA ALLISON PARK 0.7 W 40.5610 -79.9708 14.587294
## 17 US1PAAL0023 PA SOUTH PARK TOWNSHIP 0.2 NW 40.2989 -79.9970 14.648635
## 18 US1PAAL0003 PA SOUTH FAYETTE 2 SE 40.3381 -80.1159 15.392161
## 19 US1PAAL0001 PA BRIDGEVILLE 1.4 SW 40.3417 -80.1229 15.584945
## 20 US1PAAL0004 PA PENN HILLS 1.5 E 40.4759 -79.7982 16.207171
## 21 US1PAAL0039 PA ROBINSON TOWNSHIP 2.2 WNW 40.4727 -80.1656 16.406517
## 22 US1PAAL0041 PA OAKMONT 0.8 N 40.5313 -79.8345 16.682849
## 23 USC00360022 PA ACMETONIA LOCK 3 40.5361 -79.8153 18.254228
## 24 US1PAAL0006 PA MCDONALD 2.5 ENE 40.3822 -80.1871 18.323293
## 25 US1PAWS0005 PA MCMURRAY 0.2 NE 40.2831 -80.0857 18.628831
## 26 US1PAWS0008 PA VENETIA 0.7 WNW 40.2522 -80.0539 20.738838
## 27 USW00094823 PA PITTSBURGH INTL AP 40.4847 -80.2144 20.743637
## 28 US1PAAL0033 PA CARNOT-MOON 0.7 ENE 40.5224 -80.2014 21.360098
## 29 US1PAAL0030 PA CARNOT-MOON 0.9 S 40.5061 -80.2119 21.364459
## 30 USC00366111 PA MURRYSVILLE 2 SW 40.4119 -79.7244 21.730660
## 31 US1PAWT0001 PA NORTH IRWIN 2.5 WSW 40.3243 -79.7556 22.348628
## 32 USC00365918 PA MOON TOWNSHIP 40.5319 -80.2172 23.040372
## 33 USC00363343 PA GLENWILLARD DASHIELDS 40.5514 -80.2167 24.143027
## 34 USC00360355 PA BAKERSTOWN 3 WNW 40.6500 -79.9833 24.464473
## 35 US1PAWT0010 PA MURRYSVILLE 1.5 WSW 40.4317 -79.6813 25.282776
## 36 US1PAAL0012 PA SOUTH HEIGHTS 1.5 S 40.5533 -80.2379 25.760534
## 37 US1PAWS0009 PA CANONSBURG 0.3 SSW 40.2597 -80.1900 25.987112
## 38 US1PAWT0014 PA LOWER BURRELL 1.2 WNW 40.5918 -79.7333 27.543709
##
## $portugal
## [1] id name latitude longitude distance
## <0 rows> (or 0-length row.names)
##
## $quebec
## id name latitude longitude distance
## 1 CA007015417 QC NEUFCHATEL 46.8500 -71.3667 2.030466
## 2 CA007014329 QC LORETTEVILLE 46.8667 -71.3667 2.751325
## 3 CA007014298 QC LES SAULES 46.8167 -71.3167 4.105150
## 4 CA007012232 QC DUBERGER 46.8167 -71.3000 4.792654
## 5 CA007011309 QC CHARLESBOURG PARC ORLEAN 46.8667 -71.2667 5.874619
## 6 CA007011308 QC CHARLESBOURG JAR ZOO 46.9000 -71.3000 6.336821
## 7 CA007016294 QC QUEBEC/JEAN LESAGE INTL A 46.8000 -71.3833 6.462488
## 8 CA00701S001 QC QUEBEC/JEAN LESAGE INTL 46.8000 -71.3833 6.462488
## 9 CA007017BFN QC STE FOY (PIE XII) 46.7833 -71.3167 7.625681
## 10 CA007016281 QC QUEBEC 2 46.8000 -71.2667 7.874784
## 11 CA00701Q004 QC STE-FOY (U. LAVAL) 46.7833 -71.2833 8.580380
## 12 CA007012040 QC DOMAINE-DE-MAIZERETS 46.8333 -71.2167 9.560181
## 13 CA007010565 QC BEAUPORT 46.8333 -71.2000 10.808994
## 14 CA007010566 QC BEAUPORT 46.8833 -71.2000 11.269037
## 15 CA007017B65 QC STE FOY MATAPEDIA 46.7500 -71.2833 11.927703
## 16 CA007017265 QC STE FOY 46.7500 -71.2667 12.440799
## 17 CA007027700 QC ST ROMUALD 46.7667 -71.2167 13.185308
## 18 CA007018762 QC VILLENEUVE 46.8667 -71.1667 13.307164
## 19 CA007018572 QC VALCARTIER 46.9000 -71.5000 13.372471
## 20 CA007018574 QC VALCARTIER A 46.9000 -71.5000 13.372471
## 21 CA00701A9E0 QC COURVILLE DE POISSY 46.8833 -71.1667 13.685368
## 22 CA00701LF83 QC LA COURVILLOISE 46.9000 -71.1500 15.475466
## 23 CA007018577 QC VALCARTIER FES 46.9500 -71.5000 16.474749
## 24 CA007027391 QC ST JEAN CHRYSOSTOME 46.7167 -71.2167 17.545358
## 25 CA007016900 QC ST AUGUSTIN 46.7333 -71.5000 17.797745
## 26 CA007024254 QC LAUZON 46.8167 -71.1000 18.628731
## 27 CA007016932 QC STE CATHERINE 46.8500 -71.6167 21.042305
## 28 CA007047442 QC ST LAURENT ILE D'ORLEANS 46.8667 -71.0500 22.128358
## 29 CA007012240 QC DUCHESNAY 46.8667 -71.6500 23.644047
## 30 CA007020567 QC BEAUSEJOUR 46.6667 -71.1667 24.283856
## 31 CA007041330 QC CHATEAU RICHER 46.9667 -71.0333 26.668362
##
## $`rapid city`
## id name latitude longitude distance
## 1 US1SDPN0044 SD RAPID CITY 2.0 SSE 44.0487 -103.2173 1.274894
## 2 USC00396948 SD RAPID CITY WFO 44.0728 -103.2108 1.601898
## 3 US1SDPN0036 SD RAPID CITY 2.9 SE 44.0414 -103.1994 2.643470
## 4 US1SDPN0001 SD RAPID CITY 1.7 NW 44.0850 -103.2339 2.993462
## 5 US1SDPN0053 SD RAPID CITY 2.6 SW 44.0519 -103.2739 4.400391
## 6 US1SDPN0049 SD RAPID CITY 3.8 S 44.0208 -103.2418 4.694241
## 7 US1SDPN0030 SD RAPID CITY 3.3 SSW 44.0345 -103.2677 4.751190
## 8 US1SDPN0033 SD RAPID CITY 2.5 WSW 44.0644 -103.2805 4.858829
## 9 US1SDPN0043 SD RAPID CITY 3.9 S 44.0183 -103.2395 4.891807
## 10 US1SDPN0005 SD RAPID CITY 4.3 ESE 44.0472 -103.1555 5.347390
## 11 US1SDPN0027 SD RAPID CITY 2.4 WNW 44.0849 -103.2804 5.563235
## 12 US1SDPN0006 SD RAPID CITY 4.5 E 44.0661 -103.1436 6.142072
## 13 US1SDPN0057 SD RAPID CITY 4.7 E 44.0624 -103.1398 6.413886
## 14 US1SDPN0058 SD RAPID CITY 4.6 E 44.0722 -103.1402 6.518563
## 15 US1SDPN0038 SD RAPID CITY 3.5 NE 44.1134 -103.1860 6.529314
## 16 US1SDPN0002 SD RAPID CITY 3.4 WNW 44.0886 -103.2997 7.117009
## 17 US1SDPN0017 SD RAPID CITY 3.5 WNW 44.0860 -103.3015 7.123918
## 18 USC00396947 SD RAPID CITY 4NW 44.1150 -103.2828 7.909484
## 19 US1SDPN0047 SD RAPID CITY 6.3 SE 44.0061 -103.1506 8.167081
## 20 US1SDPN0004 SD RAPID CITY 5.4 SW 44.0135 -103.3008 8.273610
## 21 US1SDPN0031 SD RAPID CITY 5.8 SW 44.0057 -103.2975 8.651065
## 22 US1SDPN0016 SD RAPID CITY 6.1 SW 44.0069 -103.3094 9.270343
## 23 US1SDPN0032 SD RAPID CITY 6.3 SW 44.0028 -103.3105 9.633226
## 24 US1SDPN0029 SD RAPID CITY 6.4 WNW 44.1009 -103.3574 11.880241
## 25 US1SDPN0040 SD RAPID CITY 6.9 W 44.0866 -103.3701 12.350597
## 26 US1SDMD0019 SD BLACKHAWK 0.1 WNW 44.1520 -103.3165 12.806932
## 27 USW00024090 SD RAPID CITY RGNL AP 44.0433 -103.0536 13.427258
## 28 USW00024006 SD ELLSWORTH AFB 44.1500 -103.1000 13.854769
## 29 US1SDPN0034 SD RAPID CITY 8.1 WNW 44.1294 -103.3783 14.810927
## 30 US1SDPN0018 SD RAPID CITY 9.8 SW 43.9660 -103.3603 15.334040
## 31 US1SDMD0004 SD BLACKHAWK 2.0 NNW 44.1780 -103.3326 15.904487
## 32 US1SDPN0050 SD BOX ELDER 1.6 NE 44.1387 -103.0532 15.937019
## 33 US1SDMD0018 SD BLACKHAWK 2.3 NNW 44.1825 -103.3347 16.412472
## 34 US1SDPN0056 SD RAPID CITY 9.6 WNW 44.1274 -103.4116 17.038100
## 35 USC00394343 SD JOHNSON SIDING 44.0839 -103.4342 17.317536
## 36 USR0000SBAK SD BAKER PARK SOUTH DAKOTA 43.9792 -103.4250 18.692682
## 37 US1SDPN0015 SD HERMOSA 4.4 NNE 43.8980 -103.1560 18.727341
## 38 US1SDMD0032 SD PIEDMONT 4.6 E 44.2307 -103.2966 19.940748
## 39 USC00396427 SD PACTOLA DAM 44.0622 -103.4819 20.928415
## 40 US1SDPN0022 SD BOX ELDER 6.2 ESE 44.1000 -102.9500 22.021177
## 41 USC00394556 SD KEYSTONE 43.9039 -103.4100 23.073539
## 42 US1SDCS0018 SD HERMOSA 4.9 W 43.8494 -103.2883 24.047399
## 43 US1SDMD0001 SD BOX ELDER 8.0 NNE 44.2318 -103.0345 24.166211
## 44 US1SDCS0017 SD HERMOSA 0.3 ENE 43.8410 -103.1854 24.508709
## 45 US1SDMD0025 SD PIEDMONT 1.0 NW 44.2405 -103.4041 24.871196
## 46 USC00392888 SD FARMINGDALE 4 N 44.0333 -102.9000 25.747413
## 47 USR0000SNEM SD NEMO SOUTH DAKOTA 44.1917 -103.5097 27.370234
## 48 US1SDPN0013 SD SILVER CITY 44.0848 -103.5634 27.572203
## 49 USC00395870 SD MT RUSHMORE NATL MEM 43.8769 -103.4578 27.869317
## 50 USR0000SMRU SD MT. RUSHMORE SOUTH DAKOTA 43.8750 -103.4583 28.051427
## 51 USC00393775 SD HERMOSA 3 SSW 43.8069 -103.2131 28.148859
## 52 US1SDPN0055 SD HILL CITY 1.0 ENE 43.9382 -103.5542 29.967107
##
## $sarpsborg
## id name latitude longitude distance
## 1 NOE00109849 SARPSBORG 59.2856 11.1144 2.050694
## 2 NOE00109804 BATEROD 59.3078 11.1314 4.251261
## 3 NOE00134298 FLOTER 59.4964 11.0131 24.358920
## 4 NOE00100575 HALDEN 59.1225 11.3883 24.795447
## 5 NOE00109786 HVALER 59.0358 11.0517 27.201683
## 6 NOE00109876 MOSS BRANNSTASJON 59.4428 10.6842 28.822789
## 7 NOE00109867 MOSS 59.4339 10.6667 29.008867
##
## $`stoke-on-trent`
## [1] id name latitude longitude distance
## <0 rows> (or 0-length row.names)
##
## $sydney
## id name latitude longitude
## 1 ASN00067072 FAIRFIELD HEIGHTS POST OFFICE -33.8667 150.9333
## 2 ASN00067019 PROSPECT RESERVOIR -33.8193 150.9127
## 3 ASN00067017 GREYSTANES (BATHURST STREET) -33.8136 150.9392
## 4 ASN00067070 MERRYLANDS (WELSFORD STREET) -33.8269 150.9767
## 5 ASN00067008 GUILDFORD -33.8667 150.9833
## 6 ASN00067091 CABRAMATTA -33.9000 150.9167
## 7 ASN00067006 FAIRFIELD MWSDB -33.9000 150.9500
## 8 ASN00067032 WESTMEAD AUSTRAL AVENUE -33.8144 150.9833
## 9 ASN00067114 ABBOTSBURY (FAIRFIELD CITY FAR -33.8667 150.8611
## 10 ASN00066169 VILLAWOOD ARCHIVES -33.8333 151.0000
## 11 ASN00067119 HORSLEY PARK EQUESTRIAN CENTRE -33.8511 150.8567
## 12 ASN00067110 SEVEN HILLS (RADIO FM 103.2) -33.7858 150.9236
## 13 ASN00067013 LIVERPOOL -33.9167 150.9333
## 14 ASN00066121 CHESTER HILL -33.8833 151.0000
## 15 ASN00066025 WARWICK FARM (LIVERPOOL TREATM -33.9194 150.9383
## 16 ASN00066085 GRANVILLE RSL BOWLING CLUB -33.8361 151.0128
## 17 ASN00067035 LIVERPOOL(WHITLAM CENTRE) -33.9272 150.9128
## 18 ASN00067026 SEVEN HILLS (COLLINS ST) -33.7704 150.9318
## 19 ASN00067020 LIVERPOOL (MICHAEL WENDEN CENT -33.9214 150.8861
## 20 ASN00067064 CECIL PARK ANDERSONS RES.FARM -33.8667 150.8333
## 21 ASN00066137 BANKSTOWN AIRPORT AWS -33.9181 150.9864
## 22 ASN00067080 WINSTON HILLS LANHAMS ROAD -33.7774 150.9779
## 23 ASN00067000 EASTERN CREEK (WONDERLAND) -33.8022 150.8479
## 24 ASN00067069 GREENVALLEY (MILLER) -33.9167 150.8667
## 25 ASN00066055 LIDCOMBE (CARNARVON GOLF CLUB) -33.8667 151.0333
## 26 ASN00066134 GRANVILLE SHELL REFINERY -33.8322 151.0340
## 27 ASN00067059 BLACKTOWN -33.7694 150.8856
## 28 ASN00066168 MILPERRA BRIDGE (GEORGES RIVER -33.9289 150.9831
## 29 ASN00067042 KINGS LANGLEY (SOLANDER RD) -33.7610 150.9498
## 30 ASN00067111 NORTH PARRAMATTA (BURNSIDE HOM -33.7931 151.0167
## 31 ASN00066174 DUNDAS -33.8167 151.0333
## 32 ASN00067109 BAULKHAM HILLS EUCALYPTUS CT -33.7678 150.9814
## 33 ASN00066124 PARRAMATTA NORTH (MASONS DRIVE -33.7917 151.0181
## 34 ASN00067016 MINCHINBURY -33.8000 150.8333
## 35 ASN00066050 POTTS HILL RESERVOIR -33.8933 151.0346
## 36 ASN00066003 BANKSTOWN (CONDELL PARK) -33.9167 151.0167
## 37 ASN00067097 PRESTONS BERNERA ROAD -33.9333 150.8667
## 38 ASN00066154 HOLDSWORTHY AIR CAVALRY -33.9500 150.9500
## 39 ASN00066171 MOOREBARK N.B.GOLF CLUB -33.9500 150.9500
## 40 ASN00067120 HOXTON PARK (RANIERI PLACE) -33.9319 150.8553
## 41 ASN00066164 ROOKWOOD (HAWTHORNE AVE) -33.8771 151.0577
## 42 ASN00067112 NORTH ROCKS (MUIRFIELD GOLF CL -33.7672 151.0186
## 43 ASN00066195 SYDNEY OLYMPIC PARK (SYDNEY OL -33.8521 151.0646
## 44 ASN00066070 STRATHFIELD GOLF CLUB -33.8805 151.0631
## 45 ASN00066135 SILVERWATER (RANAD NEWINGTON) -33.8333 151.0667
## 46 ASN00066054 REVESBY (PATEN STREET) -33.9474 151.0065
## 47 ASN00067009 GLENFIELD (MACQUARIE) -33.9667 150.9000
## 48 ASN00067076 QUAKERS HILL TREATMENT WORKS -33.7366 150.8758
## 49 ASN00067092 QUAKERS HILL DOUGLAS RD. -33.7333 150.8833
## 50 ASN00067036 AUSTRAL EIGHTH AVE -33.9333 150.8167
## 51 ASN00066185 CARLINGFORD (BARELLAN AV) -33.7801 151.0587
## 52 ASN00066082 CONCORD WEST PLASTER MILLS -33.8333 151.0833
## 53 ASN00066191 GLENFIELD (HARROW ROAD) -33.9770 150.9042
## 54 ASN00067117 HOLSWORTHY CONTROL RANGE -33.9795 150.9254
## 55 ASN00067118 OAKHURST (LAWTON PLACE) -33.7431 150.8356
## 56 ASN00067102 ST CLAIR (JUBA CLOSE) -33.8044 150.7778
## 57 ASN00067100 CASTLE HILL (KATHLEEN AVE) -33.7260 150.9944
## 58 ASN00067083 MOUNT DRUITT FRANCIS STREET -33.7667 150.8000
## 59 ASN00067089 WEST PENNANT HILLS (CUMBERLAND -33.7459 151.0402
## 60 ASN00067003 COLYTON (CARPENTER ST) -33.7770 150.7877
## 61 ASN00067098 WEST PENNANT HILLS (ORATAVA A -33.7487 151.0449
## 62 ASN00066013 CONCORD GOLF CLUB -33.8523 151.0985
## 63 ASN00066087 EASTWOOD BOWLING CLUB -33.7933 151.0862
## 64 ASN00066064 CONCORD WALKER HOSPITAL -33.8333 151.1000
## 65 ASN00066076 WILEY PARK (ROSELANDS) -33.9367 151.0700
## 66 ASN00066057 RYDE PUMPING STATION -33.8167 151.1000
## 67 ASN00067061 ROSSMORE (SOUTH CREEK) -33.9353 150.7819
## 68 ASN00066048 CONCORD (BRAYS RD) -33.8483 151.1105
## 69 ASN00066020 EPPING CHESTER STREET -33.7687 151.0852
## 70 ASN00066071 GLADESVILLE CHAMPION RD -33.8381 151.1161
## 71 ASN00066017 FIVE DOCK (BARNWELL PARK GOLF -33.8681 151.1190
## 72 ASN00067024 ST MARYS BOWLING CLUB -33.7667 150.7667
## 73 ASN00067037 SCHOFIELDS BOUNDARY RD -33.6947 150.8868
## 74 ASN00066194 CANTERBURY RACECOURSE AWS -33.9057 151.1134
## 75 ASN00066144 PEAKHURST FOREST ROAD -33.9667 151.0667
## 76 ASN00066148 PEAKHURST GOLF CLUB -33.9700 151.0638
## 77 ASN00066034 ABBOTSFORD (BLACKWALL POINT RD -33.8507 151.1295
## 78 ASN00067116 WILLMOT (RESOLUTION AVE) -33.7231 150.7997
## 79 ASN00066156 MACQUARIE PARK (WILLANDRA VILL -33.7791 151.1121
## 80 ASN00067060 LONG POINT -34.0167 150.9000
## 81 ASN00067068 BADGERYS CREEK MCMASTERS F.STN -33.8683 150.7278
## 82 ASN00066047 PENNANT HILLS (YARRARA ROAD) -33.7324 151.0767
## 83 ASN00067107 VARROVILLE (ST JAMES ROAD) -33.9928 150.8178
## 84 ASN00066190 INGLEBURN (SACKVILLE STREET) -34.0117 150.8647
## 85 ASN00066049 PENSHURST -33.9667 151.0833
## 86 ASN00066000 ASHFIELD BOWLING CLUB -33.8850 151.1344
## 87 ASN00067108 BADGERYS CREEK AWS -33.8969 150.7281
## 88 ASN00066181 OATLEY (WORONORA PARADE) -33.9766 151.0766
## 89 ASN00066004 BEXLEY BOWLING CLUB -33.9430 151.1098
## 90 ASN00066018 EARLWOOD BOWLING CLUB -33.9333 151.1167
## 91 ASN00066081 NORTH RYDE STROUD STREET -33.8000 151.1333
## 92 ASN00066193 MARSFIELD (MACQUARIE UNIVERSIT -33.7654 151.1176
## 93 ASN00067025 ST MARYS MWSDB -33.7333 150.7667
## 94 ASN00067086 DURAL (OLD NORTHERN ROAD) -33.6867 151.0250
## 95 ASN00067071 THORNLEIGH BRIDGEVIEW CRESCENT -33.7167 151.0833
## 96 ASN00066036 MARRICKVILLE GOLF CLUB -33.9186 151.1402
## 97 ASN00066092 DURAL -33.6833 151.0333
## 98 ASN00066131 RIVERVIEW OBSERVATORY -33.8258 151.1556
## 99 ASN00067104 BOX HILL (HYNDS ROAD) -33.6617 150.9000
## 100 ASN00067101 BOX HILL JUNCTION ROAD -33.6649 150.8780
## 101 ASN00066189 WEST PYMBLE (WYUNA ROAD) -33.7693 151.1380
## 102 ASN00066028 HORNSBY (PRETORIA PARADE) -33.7083 151.0839
## 103 ASN00067084 ORCHARD HILLS TREATMENT WORKS -33.8020 150.7069
## 104 ASN00068043 MINTO SURREY STREET -34.0283 150.8433
## 105 ASN00066204 OYSTER BAY (GREEN POINT ROAD) -34.0009 151.0738
## 106 ASN00066032 LINDFIELD WEST -33.7822 151.1486
## 107 ASN00066069 HURSTVILLE GROVE (WAITARA PARA -33.9833 151.1000
## 108 ASN00066074 ROCKDALE BOWLING CLUB -33.9500 151.1333
## 109 ASN00066158 TURRAMURRA (KISSING POINT ROAD -33.7366 151.1271
## 110 ASN00066063 WAHROONGA RESERVOIR -33.7206 151.1128
## 111 ASN00066120 GORDON GOLF CLUB -33.7617 151.1462
## 112 ASN00068220 MINTO (ALDERNEY STREET) -34.0411 150.8458
## 113 ASN00066023 LUCAS HEIGHTS AWS -34.0514 150.9800
## 114 ASN00066078 LUCAS HEIGHTS (ANSTO) -34.0517 150.9800
## 115 ASN00066027 HORNSBY MWSDB -33.7000 151.1000
## 116 ASN00066159 HORNSBY (MOUNT WILGA) -33.6939 151.0925
## 117 ASN00067015 BRINGELLY (MARYLAND) -33.9696 150.7250
## 118 ASN00068160 CAMPBELLTOWN (KENTLYN (GEORGES -34.0542 150.8772
## 119 ASN00067106 BERKSHIRE PARK FIRST RD -33.6708 150.7972
## 120 ASN00068250 CAMDEN VALLEY GOLF RESORT -34.0128 150.7675
## 121 ASN00066060 SUTHERLAND MWSDB -34.0333 151.0667
## 122 ASN00066157 PYMBLE (CANISIUS COLLEGE) -33.7371 151.1521
## 123 ASN00066058 SANS SOUCI (PUBLIC SCHOOL) -33.9942 151.1292
## 124 ASN00066205 WAHROONGA (BOUNDARY ROAD) -33.7073 151.1315
## 125 ASN00066056 ROSEVILLE BOWLING CLUB -33.7833 151.1833
## 126 ASN00067022 GALSTON (ROWLAND VILLAGE) -33.6550 151.0553
## 127 ASN00066090 ENGADINE (CAMBRAI AVE) -34.0608 151.0125
## 128 ASN00066075 WAVERTON BOWLING CLUB -33.8411 151.1967
## 129 ASN00066010 CHATSWOOD COUNCIL DEPOT -33.8014 151.1917
## 130 ASN00068231 RUSE (DENISON STREET) -34.0630 150.8489
## 131 ASN00066114 NORTH TURRAMURRA (DRYDEN RD) -33.7179 151.1470
## 132 ASN00066037 SYDNEY AIRPORT AMO -33.9465 151.1731
## 133 ASN00066033 ALEXANDRIA (HENDERSON ROAD) -33.8992 151.1933
## 134 ASN00066061 SYDNEY NTH BOWLING CLUB -33.8333 151.2000
## 135 ASN00066067 WOLLSTONECRAFT -33.8333 151.2000
## 136 ASN00066062 SYDNEY (OBSERVATORY HILL) -33.8607 151.2050
## 137 ASN00066011 CHATSWOOD BOWLING CLUB -33.8000 151.2000
## 138 ASN00066012 CHATSWOOD WATER SUPPLY -33.8000 151.2000
## 139 ASN00067018 PENRITH LADBURY AVENUE -33.7542 150.6783
## 140 ASN00067073 MARALYA BOUNDARY ROAD -33.6200 150.8967
## 141 ASN00067115 GLENMORE PARK (CARTWRIGHT CL) -33.7826 150.6619
## 142 ASN00066021 ALEXANDRIA (ERSKINEVILLE) -33.9167 151.2000
## 143 ASN00066040 MIRANDA (BLACKWOOD ST) -34.0405 151.0982
## 144 ASN00068015 CAMPBELLTOWN 2 MWSDB -34.0667 150.8167
## 145 ASN00066117 TURRAMURRA NORTH -33.7167 151.1667
## 146 ASN00066006 SYDNEY BOTANIC GARDENS -33.8662 151.2160
## 147 ASN00066080 CASTLE COVE (ROSEBRIDGE AVE) -33.7809 151.2044
## 148 ASN00066001 AUDLEY NATIONAL PARK BOTTOM ST -34.0667 151.0500
## 149 ASN00066176 AUDLEY (ROYAL NATIONAL PARK) -34.0658 151.0567
## 150 ASN00066139 PADDINGTON -33.8833 151.2167
## 151 ASN00067029 WALLACIA POST OFFICE -33.8637 150.6410
## 152 ASN00066167 NORTHBRIDGE BOWLING CLUB -33.8057 151.2145
## 153 ASN00066130 NORTHBRIDGE (SAILORS BAY) -33.8167 151.2167
## 154 ASN00067004 EMU PLAINS -33.7500 150.6667
## 155 ASN00068037 KENNY HILL -34.0500 150.7667
## 156 ASN00066166 CREMORNE GRASMERE ROAD -33.8283 151.2217
## 157 ASN00067087 GLENORIE (GATELEIGH PARK) -33.6167 151.0167
## 158 ASN00067067 EMU PLAINS -33.7600 150.6567
## 159 ASN00066206 ST IVES (RICHMOND AVENUE) -33.7096 151.1730
## 160 ASN00067113 PENRITH LAKES AWS -33.7195 150.6783
## 161 ASN00068257 CAMPBELLTOWN (MOUNT ANNAN) -34.0615 150.7735
## 162 ASN00068230 CAMPBELLTOWN U.W.S. -34.0697 150.7878
## 163 ASN00068081 CAMPBELLTOWN SWIMMING CENTRE -34.0833 150.8167
## 164 ASN00066007 BOTANY NO.1 DAM -33.9333 151.2167
## 165 ASN00066073 RANDWICK RACECOURSE -33.9105 151.2276
## 166 ASN00068254 MOUNT ANNAN BOTANIC GARDEN -34.0673 150.7678
## 167 ASN00066009 CENTENNIAL PARK ROUND HOUSE -33.8975 151.2328
## 168 ASN00067063 COBBITTY (CUTHILL) -33.9833 150.6667
## 169 ASN00066160 CENTENNIAL PARK -33.8959 151.2341
## 170 ASN00067031 WINDSOR BOWLING CLUB -33.6100 150.8151
## 171 ASN00066042 MOSMAN (BAPAUME ROAD) -33.8194 151.2428
## 172 ASN00066188 BELROSE (EVELYN PLACE) -33.7402 151.2173
## 173 ASN00068227 AMBARVALE CLENNAM AVE -34.0867 150.7917
## 174 ASN00066145 SEAFORTH CASTLE CIRCUIT -33.7894 151.2392
## 175 ASN00068207 COBBITY (ROSENEATH) -34.0167 150.6833
## 176 ASN00066184 MOSMAN COUNCIL -33.8340 151.2484
## 177 ASN00066052 RANDWICK BOWLING CLUB -33.9096 151.2419
## 178 ASN00066086 CRONULLA STP -34.0313 151.1642
## 179 ASN00066041 MOSMAN WATER SUPPLY -33.8333 151.2500
## 180 ASN00066136 CARINGBAH (LILLI PILLI) -34.0667 151.1167
## 181 ASN00068070 WORONORA DAM -34.1167 150.9333
## 182 ASN00068011 CAMDEN BOWLING CLUB -34.0500 150.7167
## 183 ASN00067010 GLENORIE (OLD NORTHERN RD) -33.5908 151.0094
## distance
## 1 1.881790
## 2 3.769152
## 3 4.135739
## 4 5.020100
## 5 5.260294
## 6 5.693717
## 7 5.858326
## 8 6.317220
## 9 6.627566
## 10 6.726345
## 11 6.770114
## 12 7.163157
## 13 7.422955
## 14 7.448599
## 15 7.754869
## 16 7.801609
## 17 8.729832
## 18 8.852678
## 19 8.913714
## 20 9.120153
## 21 9.189477
## 22 9.206123
## 23 9.260935
## 24 9.442011
## 25 9.717611
## 26 9.806921
## 27 9.856503
## 28 10.049571
## 29 10.064021
## 30 10.206744
## 31 10.234628
## 32 10.300292
## 33 10.404864
## 34 10.521499
## 35 10.790772
## 36 10.911531
## 37 10.951350
## 38 11.271660
## 39 11.271660
## 40 11.422627
## 41 12.169844
## 42 12.319787
## 43 12.431978
## 44 12.748603
## 45 12.760934
## 46 12.928587
## 47 13.268493
## 48 13.567795
## 49 13.675242
## 50 13.970001
## 51 14.205035
## 52 14.279472
## 53 14.321038
## 54 14.405998
## 55 14.744007
## 56 14.945410
## 57 15.017779
## 58 15.167677
## 59 15.416884
## 60 15.450666
## 61 15.478987
## 62 15.562401
## 63 15.746686
## 64 15.809996
## 65 16.122048
## 66 16.132769
## 67 16.638119
## 68 16.669913
## 69 16.950934
## 70 17.237878
## 71 17.567447
## 72 17.704003
## 73 17.724215
## 74 18.028240
## 75 18.097853
## 76 18.179759
## 77 18.423361
## 78 18.550317
## 79 18.579019
## 80 18.741684
## 81 18.781169
## 82 18.835558
## 83 18.955452
## 84 18.962689
## 85 19.197135
## 86 19.269018
## 87 19.355573
## 88 19.523783
## 89 19.553323
## 90 19.564367
## 91 19.585383
## 92 19.721134
## 93 19.902615
## 94 20.170027
## 95 20.504349
## 96 20.849104
## 97 20.851122
## 98 21.009524
## 99 21.120894
## 100 21.136126
## 101 21.209116
## 102 21.226965
## 103 21.288391
## 104 21.378555
## 105 21.391107
## 106 21.556418
## 107 21.581812
## 108 21.810557
## 109 22.152621
## 110 22.190861
## 111 22.258367
## 112 22.624322
## 113 22.864610
## 114 22.897282
## 115 22.914815
## 116 22.954030
## 117 23.124628
## 118 23.222412
## 119 23.404366
## 120 23.504569
## 121 23.967617
## 122 24.058869
## 123 24.391064
## 124 24.466552
## 125 24.547893
## 126 24.583503
## 127 24.644041
## 128 24.650159
## 129 24.770808
## 130 24.837610
## 131 24.858762
## 132 24.870760
## 133 24.916043
## 134 25.005231
## 135 25.005231
## 136 25.421751
## 137 25.553205
## 138 25.553205
## 139 25.580355
## 140 25.759544
## 141 25.877100
## 142 26.004107
## 143 26.257057
## 144 26.264221
## 145 26.424281
## 146 26.470162
## 147 26.489163
## 148 26.516139
## 149 26.689966
## 150 26.728506
## 151 26.729649
## 152 26.737272
## 153 26.738727
## 154 26.749858
## 155 26.859908
## 156 27.048941
## 157 27.152433
## 158 27.162530
## 159 27.351861
## 160 27.416523
## 161 27.594122
## 162 27.727260
## 163 27.966933
## 164 28.037230
## 165 28.284461
## 166 28.418732
## 167 28.449544
## 168 28.460442
## 169 28.535385
## 170 28.724326
## 171 29.091049
## 172 29.221247
## 173 29.247143
## 174 29.347935
## 175 29.352997
## 176 29.459874
## 177 29.545881
## 178 29.549535
## 179 29.612274
## 180 29.616192
## 181 29.657248
## 182 29.692767
## 183 29.742533Not all the locations have stations nearby. Therefore, I will omit them from the weather data evaluation using the following code.
has_stations <- sapply(stations, function(x) nrow(x) > 0)
outbreak_loc_true <- outbreak_loc %>%
arrange(id) %>%
filter(has_stations)
outbreak_loc_true$stations <- NULL
outbreak_loc_true## id file_id latitude longitude year_min date_min
## 1 bovenkarspel bovenkarspel 52.70 5.24 1989 1989-01-01
## 2 bronx bronx 40.82 -73.92 2005 2005-01-01
## 3 columbus columbus 39.98 -82.99 2003 2003-01-01
## 4 genesee1 genesee1 43.09 -83.63 2004 2004-01-01
## 5 genesee2 genesee2 43.09 -83.63 2005 2005-01-01
## 6 london london 51.52 -0.10 1979 1979-01-01
## 7 melbourne melbourne -37.86 145.07 1990 1990-01-01
## 8 miyazaki miyazaki 31.89 131.34 1992 1992-01-01
## 9 murcia murcia 37.98 -1.12 1991 1991-01-01
## 10 pamplona pamplona 42.81 -1.65 1996 1996-01-01
## 11 pittsburgh pittsburgh 40.43 -79.98 2002 2002-01-01
## 12 quebec quebec 46.85 -71.34 2002 2002-01-01
## 13 rapid city rapid_city 44.06 -103.22 1995 1995-01-01
## 14 sarpsborg sarpsborg 59.28 11.08 1995 1995-01-01
## 15 sydney sydney -33.85 150.93 2006 2006-01-01
## year_max date_max onset before_onset
## 1 1999 1999-12-31 1999-02-25 1999-02-11
## 2 2015 2015-12-31 2015-07-12 2015-06-28
## 3 2013 2013-12-31 2013-07-09 2013-06-25
## 4 2014 2014-12-31 2014-06-06 2014-05-23
## 5 2015 2015-12-31 2015-05-04 2015-04-20
## 6 1989 1989-12-31 1989-01-01 1988-12-18
## 7 2000 2000-12-31 2000-04-17 2000-04-03
## 8 2002 2002-12-31 2002-07-18 2002-07-04
## 9 2001 2001-12-31 2001-06-26 2001-06-12
## 10 2006 2006-12-31 2006-06-01 2006-05-18
## 11 2012 2012-12-31 2012-08-26 2012-08-12
## 12 2012 2012-12-31 2012-07-18 2012-07-04
## 13 2005 2005-12-31 2005-05-26 2005-05-12
## 14 2005 2005-12-31 2005-05-12 2005-04-28
## 15 2016 2016-12-31 2016-04-25 2016-04-11Using the countyweather codes I can gather the data for each station in a loop. The code gathers the weather data for each stations and averages them. Then I saved all the data as rds. files because they take a long time to gather. The data is saved in a folder I created called “weather_files/”
for(i in which(has_stations))
{
meteo_df <- meteo_pull_monitors(monitors = stations[[i]]$id,
keep_flags = FALSE,
date_min = outbreak_loc$date_min[i],
date_max = outbreak_loc$date_max[i],
var = c("prcp","snow","snwd","tmax","tmin","tavg"))
coverage_df <- rnoaa::meteo_coverage(meteo_df, verbose = FALSE)
filtered <- countyweather:::filter_coverage(coverage_df, 0.90)
good_monitors <- unique(filtered$id)
filtered_data <- dplyr::filter(meteo_df, id %in% good_monitors)
averaged <- countyweather:::ave_weather(filtered_data)
file_name <- paste0("weather_files/", outbreak_loc$file_id[i], ".rds")
saveRDS(averaged, file_name)
#readRDS(file_name)
}Now that all of the data is gathered and averaged I can plot the data. The loop will go through the files in order which is in alphabetical order. Therefore I must order my outbreak data frame into alphabetical order too. I will rename this data frame as df_stations for plotting.
## id file_id latitude longitude year_min date_min
## 1 bovenkarspel bovenkarspel 52.70 5.24 1989 1989-01-01
## 2 bronx bronx 40.82 -73.92 2005 2005-01-01
## 3 columbus columbus 39.98 -82.99 2003 2003-01-01
## 4 genesee1 genesee1 43.09 -83.63 2004 2004-01-01
## 5 genesee2 genesee2 43.09 -83.63 2005 2005-01-01
## 6 london london 51.52 -0.10 1979 1979-01-01
## 7 melbourne melbourne -37.86 145.07 1990 1990-01-01
## 8 miyazaki miyazaki 31.89 131.34 1992 1992-01-01
## 9 murcia murcia 37.98 -1.12 1991 1991-01-01
## 10 pamplona pamplona 42.81 -1.65 1996 1996-01-01
## 11 pittsburgh pittsburgh 40.43 -79.98 2002 2002-01-01
## 12 quebec quebec 46.85 -71.34 2002 2002-01-01
## 13 rapid city rapid_city 44.06 -103.22 1995 1995-01-01
## 14 sarpsborg sarpsborg 59.28 11.08 1995 1995-01-01
## 15 sydney sydney -33.85 150.93 2006 2006-01-01
## year_max date_max onset before_onset
## 1 1999 1999-12-31 1999-02-25 1999-02-11
## 2 2015 2015-12-31 2015-07-12 2015-06-28
## 3 2013 2013-12-31 2013-07-09 2013-06-25
## 4 2014 2014-12-31 2014-06-06 2014-05-23
## 5 2015 2015-12-31 2015-05-04 2015-04-20
## 6 1989 1989-12-31 1989-01-01 1988-12-18
## 7 2000 2000-12-31 2000-04-17 2000-04-03
## 8 2002 2002-12-31 2002-07-18 2002-07-04
## 9 2001 2001-12-31 2001-06-26 2001-06-12
## 10 2006 2006-12-31 2006-06-01 2006-05-18
## 11 2012 2012-12-31 2012-08-26 2012-08-12
## 12 2012 2012-12-31 2012-07-18 2012-07-04
## 13 2005 2005-12-31 2005-05-26 2005-05-12
## 14 2005 2005-12-31 2005-05-12 2005-04-28
## 15 2016 2016-12-31 2016-04-25 2016-04-11PLOT 1
Outbreak Distribution
This plot is divided by outbreaks in the northern and southern hemisphere. This allows us to see when the outbreaks generally occur in the year.
PLOT 2
10 years for all data
These plots allow for a quick glance into all the weather variables for each location.
## Warning: Removed 4 rows containing missing values (geom_path).
## Warning: Removed 19 rows containing missing values (geom_path).
## Warning: Removed 27 rows containing missing values (geom_path).
PLOT 3
TMAX and TMIN
## Warning: Removed 1098 rows containing missing values (geom_path).
## Warning: Removed 7 rows containing missing values (geom_path).
## Warning: Removed 3 rows containing missing values (geom_path).
## Warning: Removed 55 rows containing missing values (geom_path).
PLOT 4
Precipitation
I also made a loop to plot graphs and histograms of the data with lines indicating each day before the start of the outbreak for a total of 14 days.A plot of the percentiles is also included.
## Warning: Removed 24 rows containing missing values (geom_path).
## Warning: Removed 27 rows containing missing values (geom_path).